Temperature contributes to host specialization of coffee wilt disease (Fusarium xylarioides) on arabica and robusta coffee crops

Coffee wilt disease, caused by the fungus Fusarium xylarioides, is a vascular wilt disease that has affected coffee production in sub-Saharan Africa over the past century. Today, the disease has two host-specific populations specialising on arabica and robusta coffee crops, which grow at high and low altitude, respectively. Here we test whether adaptation to different temperatures contributes to specialisation of the fungi on each crop. Firstly, climate models show that the severity of the arabica and robusta populations of coffee wilt disease correlates with temperature. The robusta population shows higher peak severity than the arabica population overall, but the latter has greater cold tolerance. Secondly, growth assays of thermal performance of fungal strains in vitro show that, while robusta strains grow faster than arabicas at intermediate temperatures, the arabica strains have higher sporulation and spore germination rates at temperatures below 15ºC. The match between environmental patterns of severity in nature with thermal performance of fungal cultures in the laboratory supports a role for temperature adaptation in specialisation on arabica and robusta coffee. Extrapolating our temperature-models to future climate change predicts that disease severity could decline on average due to increased temperature but could increase in some coffee-growing regions.

. (a) A topographic map for elevations around the Ethiopian Highlands and the Congo Basin, with colour of the map indicating elevation. Red circles represent the 51 villages with recorded incidences of the arabica population of CWD, and blue triangles represent the 60 villages with recorded incidences of the robusta population by Oduor et al. 17 . Black numbers and lines show the localities of the four isolates used for growth assays. Isolates specialised to arabica coffee (Coffea arabica): 1 = arabica isolate 389,563, 2 = arabica isolate 393,593. Isolates specialised on 'robusta' coffee (C. canephora): 3 = robusta isolate 392,263, 4 = robusta isolate 392,278. Isolate numbers refer to IMI (International Mycological Institute) numbers in the CABI culture collection. (b) A boxplot for local elevation of villages with incidences of CWD surveyed by Oduor et al. 17 , illustrating higher elevation for villages with the arabica population of CWD.

Scientific Reports
| (2023) 13:9327 | https://doi.org/10.1038/s41598-023-36474-w www.nature.com/scientificreports/ in localities of arabica strains, which can also limit growth and transmission of the fungus, and perhaps explain in part the lower severity 20 . The effect of temperature on F. xylarioides life cycles has not been investigated previously. Studies of other Fusarium species on different crops have shown temperature has a significant impact on growth and reproduction of the fungus, however. Cruz et al. observed faster radial growth of 14 F. oxysporum isolates from soybean (FO36-FO49) on agar plates between 25 and 30 °C than at lower temperatures 21 , and this finding matched disease severity in soybean seedlings inoculated at different temperatures. Similarly, Manstretta & Rossi found that spore maturation of F. graminearum strains was restricted to a temperature of between 20 and 25 °C 22 . Past field observations have also shown CWD has higher incidence rates in warmer years 14 , which coincides with observations in other Fusarium species. Temperature is predicted to increase by 2 to 6 °C in most coffee-growing regions of Africa by 2080 23 . This could potentially facilitate the spread of the disease into previously CWD-free regions 1 .
We hypothesize that the two host-specific populations of CWD have different responses to temperature, given the distinct geographical distributions of the two 16 . Specifically, we predict that they should differ in optimal temperature (T opt , reflecting adaptation to average temperature of local climates) and critical thermal minimum (CT min , reflecting extreme low temperatures of local climates) 24 . If this is true, management strategies of arabica and robusta coffee might need to be assessed separately and adapted to maximise their efficacy against different CWD infections. Thus, to adapt cultivation and management strategies for coffee production, and cope with potential damage caused by climate change, it is necessary to understand the impacts of temperature on the growth and reproduction of strains from each F. xylarioides host-specific population.
Here, we investigate the effects of temperature on the distribution, growth and reproduction of F. xylarioides using a combination of environmental modelling and growth assays in the laboratory. Firstly, we test whether temperature predicts the distribution and severity of F. xylarioides infection using geographical records on severity, measured as the proportions of trees infected per farm, as a measure of disease intensity and of coffee plant loss 17 , and whether severity is predicted to increase with projected future climate change. Secondly, we quantify thermal performance curves of F. xylarioides isolates from arabica and robusta populations in vitro in the laboratory. Specifically, we test the hypothesis that the arabica population of F. xylarioides grows and reproduces better than the robusta population at lower temperatures, as it is adapted to the colder montane climate in the Ethiopian Highlands. We use mycelial growth, spore production and germination in vitro under a range of temperatures in the laboratory as an assay of thermal performance of growth and reproduction, of the fungus. We could not practically grow coffee plants at this range of temperatures in the laboratory to perform the assays in vivo, but we test whether evidence for thermal dependence of growth and sporulation in vitro reflects climatic correlates of distribution and severity in the field 12 . Together, the results help explain past differences in severity as well as inform predictions for how future climate change might affect the distribution and severity of disease caused by the two host-specific populations.

Materials and methods
Climatic modelling. Severity data. Severity data were extracted from the CWD survey report in Ethiopia, Uganda, and Tanzania as the proportion of trees infected by F. xylarioides per farm 17 . The mean latitude and longitude of regions/villages were extracted using ElevationMap 25 . Incidences in Ethiopia were assigned as the arabica population (n = 51), and those in Uganda and Tanzania as the robusta population (n = 60), as stated by Oduor et al. 17 (Fig. 1a). The location and severity records of CWD were loaded using R packages 'readr' and 'dismo' 26,27 .
Climate data. Bioclimatic variables, which are measures of global temperature and precipitation patterns from 1970 to 2000, were downloaded from the WorldClim database at a spatial resolution of 10 min 23 . As samples of CWD severity were taken during 2002 17 , we extracted data for 2002 from the WorldClim database and calculated bioclimatic variables for this year 28 to improve accuracy of prediction. CMIP6-projected bioclimatic variables from 2060 to 2080 were also downloaded at a spatial resolution of 10 min 29 . We chose the shared socioeconomic pathway of SSP2-4.5 using predictions of the HadGEM3-GC3.1 model for the CMIP6 data, which represent the intermediate scenario of projected future climate change 29 . The data of global bioclimatic variables were downloaded and transformed for data analysis using R packages 'dismo' and 'raster' 27,30 .
Statistical modelling. To assess the effect of temperature on severity of CWD, we extracted bioclimatic temperature variables of villages with incidences of the disease in 2002, using the latitude and longitude data in Oduor et al. 17 , with R package 'raster' 30 . We expected the transmission of CWD to have a non-linear thermal response, like other fungal plant diseases 31 . Thus, we fitted the severity of arabica and robusta population of F. xylarioides as response variables using the Sharpe-Schoolfield equation 32,33 , which is derived from the thermal sensitivity of biochemical reactions and has been widely used for thermal performance curves 34 . This model fits the functional response of trait R (in this case severity) to temperature as where R ref is the trait value at a standard temperature T ref , E is the activation energy (eV), which describes how quickly performance drops at temperatures below the optimum, E D is the high temperature de-activation energy (eV), which describes how quickly performance drops at temperatures above the optimum, T h is the temperature www.nature.com/scientificreports/ where half of the enzyme units underlying the trait are assumed to be inactive, T is the temperature in ºC and k is the Boltzmann constant (= 8.62 × 10 -5 eV/K). We used bio1 (annual mean temperature), bio5 (highest temperature of hottest month), and bio6 (lowest temperature of coldest month) as explanatory variables in the models to investigate the response of CWD transmission to average and extreme temperatures respectively. The model was fitted using non-linear least squares (NLS) with the 'sharpeschoolhigh_1981' function R packages 'rTPC' , 'nls.multstart' and 'broom' 32,33 and using the 'get_start_vals' function to obtain sensible starting values. The fitted models estimated the optimum temperature (T opt ), severity at optimum temperature (R max ), and critical thermal minimum (CT min ) and maximum (CT max ) of the populations.
Predicting potential severity in other regions and into the future. We used the fitted models to predict potential severity of arabica and robusta populations across all African coffee-producing countries using WorldClim  bioclimatic data, to reflect the period leading up to the original survey by Oduor et al. 17 . We predicted potential future severity under projected climate changes for the region using CMIP6 (2060-2080) bioclimatic data by predicting values for the model parameters but fed with future temperatures instead of past ones. As temperature will not be the only variable affecting the severity of CWD, we also reported models fitting additional bioclim variables to incorporate precipitation and other seasonal and daily variables (Supplementary Information S1), although our focus is primarily on temperature here.
Growth assays. Isolate selection. All isolates used in the experiments were provided by the CABI-IMI culture collection. The isolates used in this study were accessed and used in full compliance with UK Regulation implementing the Nagoya Protocol on Access and Benefit Sharing https:// www. gov. uk/ guida nce/ abs of the Convention on Biodiversity. Four F. xylarioides isolates, two from the arabica and two from the robusta populations, were chosen from locations with typical montane and tropical climatic patterns respectively (Fig. 1a). Specifically, average monthly temperature and precipitation of each location was extracted from the WorldClim database 23 and cross-referenced to pick isolates from locations spanning the range of monthly temperatures and precipitations observed for that host population. Isolates were used for growth assays and spore production and germination measurements. (Table S1).
Temperature selection. Based on the current and projected range of temperatures experienced using World-Clim and CMIP6 records 23,29 , we conducted growth assays at 10, 15, 20, 25, 30, 35, and 40 °C to analyse response in growth and sporulation rates. Three replicates of each F. xylarioides isolate was grown at each temperature to measure variation in responses. Due to the limited availability of incubators, growth assays were performed in 2 batches. The first round included 15, 20, 25, 30, and 40 °C, while the second covered 10, 15, 20, 25, 30 and 35 °C to ensure thermal performance in growth and reproduction of isolates were captured under a range of temperatures. We conducted the first batch of growth assays from 20th to 31st January 2021, and the second batch from 7 to 18th May 2021. Note that the 40 °C and 35 °C treatments were only conducted in batch 1 and 2, respectively, 10 °C was also conducted only in batch 2, whereas all other treatments were repeated in both batches.
Growth assays. Assays were performed on 120 mm synthetic nutrient agar (SNA) plates based on the protocol of Halder 35 . A block of agar, approximately 3 mm × 3 mm × 3 mm, was transferred from the active edge of fungal cultures (stock plates); these stock plates contained mycelium of the corresponding isolate and inocula were placed at the centre of a fresh plate using a sterile disposable loop. A fine-tip marker was used to trace the outline of the agar block on the base of each 60mm Petri dish. The plates were grown at 25 °C for 48 h, after which the new colony outline was drawn in a different colour and photographed. Plates were then randomly assigned to and grown at the range of temperatures for 9 days (216 h, 3 replicates for each strain x temperature x batch combination). Photographs were taken of each plate every 24 h and analysed using Photoshop (Adobe Photoshop CC 2015.5) and ImageJ (ImageJ 1.53c) to count the pixel numbers of the fungal colony and the plate. The growth rate of the colony was then calculated using the following equations: Average daily growth rate between days four and eight was used to represent growth rates, reflecting the period of maximum growth after an initial lag and before a final slow down.
Spore counts. After the last sampling from the growth assay was completed, fungal spores were harvested from two out of three replicates per isolate per temperature, by pouring 10 mL distilled water and 1 droplet of 0.05% tween-80 (polysorbate 80) solution onto the mycelia and scraping its surface to release spores from the colony. Spore concentration was measured for each sample using a haemocytometer, with the average spore concentration in the plate calculated following BTI protocol 36 . Each of the two replicates from each isolate per temperature was sampled twice, with 1 measurement for each of 10 haemocytometer chambers. This gave a total of 20 repeated measurements of spore concentration for each replicate, which were averaged. We then calculated sporulation rate as the number of spores produced per mm 2 of mycelium by dividing the average spore concen- www.nature.com/scientificreports/ tration × 10 mL by the area of the mycelium for that sample. This measure aimed to account for differential spore production rate due to varied sizes of colonies.
Spore germination rate. This was measured by spreading a thin line of inoculum produced from the last section in a fresh SNA plate. Each combination of strain × temperature × batch had 1 plate sampled. The plates were kept at room temperature for 24 h to allow the spores to germinate and then viewed under a compound microscope (Meiji ML2000 microscope) with × 4 magnification. The total number of spores in the sight and the number of germinated spores with hyphae growing out were recorded. Counts were made from 3 different points on each plate. Germination rate was calculated as the number of germinated spores divided by the total number of spores, averaged across the 3 points measured on each plate.
Statistical analyses. We used a linear mixed-effects model (LMM), with average daily growth rate as the response variable, and factorised temperature, host-specific population and their interactions as fixed effects, and batch and strain as random effects, using R packages 'lme4' and 'lmerTest' 37,38 . By analysing temperature as factors, we could compare nonlinear responses to temperature within and between populations. Pairwise comparisons between mean growth rates at different temperatures were conducted with Tukey adjustments using R package 'emmeans' 39 . The same approach was then repeated for sporulation rates and spore germination rates as response variables in turn. A NLS model was also fitted to the three response variables in turn (growth rate, sporulation rate, germination rate) at each incubation temperature, following the Sharpe-Schoolfield model (high-temperature inactivation only) as described above for climate data. We estimated T opt , R max , CT max and CT min for each response variable in turn. Code and data for all statistical analyses is available at github.

Results
Climatic modelling of the thermal dependence of severity of Fusarium xylarioides. Response to mean annual temperature. The non-linear least squared (NLS) model fitted a flat thermal response curve for the arabica host-specific populations, with low severity across a wider range of temperatures. In contrast, the severity of the robusta host-specific population of F. xylarioides showed a strong bell-shaped response to mean annual temperatures, with higher optimum mean temperatures and maximum severities for the than the arabica form (Fig. 2) . The model predictions also showed the arabica population displays higher severity than the robusta population below approximately 20.5 °C, while the latter has higher severity above this mean temperature (Fig. 2a). The only cases with severity above 10% for arabica were in the lower half of the observed temperature range.
Response to extreme temperatures. NLS models for extreme temperatures also predicted the arabica population to have lower optimum temperature, with better tolerance to coldness and higher severity at low temperatures. The models for maximum temperature suggested the severity of the arabica population to be highest ( Despite the robusta population having much higher peak severity, the arabica population is predicted to have higher severity than the former at locations with maximum temperature below approximately 27.2 °C (Fig. 2b).
The minimum temperature models also revealed a similar difference between the populations. Severity of the arabica population was predicted to peak (3.41%, 95% CL 2.11-14.10%) at minimum temperature of 4.2 °C (95% CL 3.0 °C-16.2 °C), with CT min outside the range of minimum temperatures across sampled sides (95% CL -Inf-5.7 °C). In contrast, peak severity of the robusta population of 45.2% (95% CL 36.6-51.0%) was predicted at minimum temperature of 17.1 °C (95% CL 16.9-17.5 °C), with CT min at 14.36 °C (95% CL 12.7-14.5 °C). The arabica population is predicted to have higher severity than robusta below minimum annual temperature of 14 °C (Fig. 2c).

Spatial predictions of severity of Fusarium xylarioides. Past severity. WorldClim (1970-2000) data
were used to extrapolate the models reported above to predict past spatial variation in severity across coffeeproducing regions of Africa (between 25°S and 20°N, distribution of coffee producing countries can be found in Fig. S1) based on temperature alone. The arabica population is confined to Ethiopia, and its predicted severity is, as expected, also low across Africa under 1970-2000 climate for all three measures of temperature (mean temperature model: mean = 1.12%, SD = 0.690%; maximum temperature model: mean = 1.14%, SD = 0.996%; minimum temperature model: mean = 1.40%, SD = 0.914%). Some coffee-farming regions were shown to have higher predicted severity than other regions (over 4%), however, including the Ethiopian mountains (the current distribution of the arabica population), Burundi, Rwanda, southern Kenya, and Zimbabwe (all regarded as important arabica coffee-producing regions by ICO 40 (Fig. 3a,b,c). The mean and maximum temperature models also projected Tanzania as a suitable habitat for the arabica isolates (Fig. 3a,b), whereas the minimum temperature model included western Zambia (Fig. 3c). Although not a major producer of arabica coffee, the models also predict Angola and northern South Africa to be suitable habitats of the arabica population (with predicted severity > 4%, Fig. 3a,b,c). www.nature.com/scientificreports/ All three models for the robusta population predicted on average higher severity compared to the corresponding models of the arabica population (mean temperature model: mean = 8.89%, SD = 11.172%; maximum temperature model: mean = 23.52%, SD = 15.256%; minimum temperature model: mean = 7.86%, SD = 12.642%). The robusta coffee-growing regions predicted to have high severity (> 30%) include northern Tanzania, the Congo Basin (including Uganda, Democratic Republic of the Congo, Central African Republic, and Republic of the Congo), Cameroon, west Africa (including Côte d'Ivoire, Gabon and Guinea), and Nigeria (Fig. 4a,b,c) 14,41 . Furthermore, the maximum-temperature model predicted most regions south of 10°N to have high severity (Fig. 4b). Though Malawi, southern Ethiopia and Madagascar are not major producers of robusta coffee, these countries were also estimated to be climatically favourable for high severity of the robusta population by all three models (Fig. 4a,b,c).
Future severity. CMIP6 (2060-2080) bioclimatic data were used to extrapolate the models to predict future potential severity across coffee-producing regions of Africa based on future temperature projections, assuming correlation of severity with temperature variables remained unchanged. All three models predicted a significant near-half reduction in severity of the arabica population under projected 2060-2080 temperature changes (Welch Two Sample t-test, mean temperature model: mean change = − 0.52%, t 108,089 = − 165.41; maximum temperature model: mean change = − 0.61%, t 106,150 = − 137.36; minimum temperature model: mean change = − 0.53%, t 126,142 = − 114.72. P-value < 0.001 for all three models). However, some regions near Bale Mountains, Simien Mountains and Ethiopian Highlands were shown to experience an increase in severity of CWD by all models, despite the overall decrease in severity in this region (> 2%, Fig. 3d,e,f). In addition, arabica population severity was predicted to increase in some locations of south Kenya, Rwanda, Burundi and east DRC, by the mean and maximum temperature models (Fig. 3d,e).
The mean and maximum-temperature models also predicted a near-half reduction in severity of the robusta population under 2060-2080 climate, while the minimum-temperature model predicted a slight increase (Welch Two Sample t-test, mean temperature model: mean change = − 4.51%, t 128,589 = − 79.30; However, in all models, some regions would experience a large increase in severity (> 20%), and the major robusta-coffee-producing areas subjected to this predicted increase were east Democratic Republic of the Congo, Tanzania, and Uganda (Fig. 4d,e,f) 40 . Angola and Madagascar also produce some robusta coffee 41 , and were predicted to have increased severity (Fig. 4d,e). Rwanda and Burundi were predicted to experience increased severity, but they are currently arabica-producing regions. As context for our temperature-only models, climate models incorporating additional bioclimatic variables yielded similar average levels for net changes in severity by 2060-2080, but the exact regions projected to  www.nature.com/scientificreports/ experience increased severity were different from the temperature-only models for both host-specific populations ( Supplementary Information S1, Fig. S2).
Growth assays in the laboratory. Mycelial growth rates. The linear mixed-effect model accounted for a large proportion of variations in growth rates in replicates across two F. xylarioides isolates of each host-specific population (model adjusted R 2 = 0.777), and the fixed effects accounted for most of variation (fixed effects adjusted R 2 = 0.771). There was a significant interaction between host and temperature (F 6,102 = 4.9, P < 0.001, Fig. 5a, Table S2) and the growth rates of the robusta population significantly exceeded the arabica population at 20 °C . Contrary to our prediction that arabica isolates should grow significantly faster at low temperatures, isolates of both populations showed a similar shape of thermal response: slow growth rates at 10 °C and 15 °C, a steep increase at both 20 °C and 25 °C, and a slower decline at 30 °C and above (model outputs in Table S2). The Sharp-Schoolfield model predicted an optimum temperature for growth of arabica strains at 22.3 °C (95% CL 20.7-23.6) with maximum growth rate of 2.14 mm/day (95% CL 1.9-2.5) and CT min at 12.0 °C (95% CL 8.2-15.0, Fig. 5a). The robusta population was also predicted to grow fastest at 22.3 °C (95% CL 21.5-22.9) with growth rate of 3.72 mm/day (95% CL 3.3-4.4) and CT min at 14.5 °C (95% CL 11.5-17.4). CT max in both cases was extrapolated well outside experimental temperatures (> 50 °C) and hence not reported. The fitted model predicted that arabica isolates grow slightly faster than robusta isolates below 17 °C, but this difference was not significant in the LME model.
Spore germination rates. The linear mixed-effects model explained most of the variation in spore germination rates of F. xylarioides isolates (model adjusted R 2 = 0.725), and temperature and population accounted for most of the variation (fixed effects adjusted R 2 = 0.512). Due to lower replication for the germination assay, the interaction between host population and temperature was not significant (F 6,27 1.1, P = 0.39, Fig. 5c, Table S4) but the pattern was qualitatively similar to sporulation rate: the populations showed similar germination rates at higher temperatures, but arabica strains had higher germination rates than robusta strains at lower temperatures (Fig. 5c) Fig. 5c).

Discussion
In both climate modelling and growth assays, the two host-specific populations of F. xylarioides showed distinct responses to temperature that reflect their elevational distributions. The NLS Sharp-Schoolfield models imply the optimum temperatures (T opt ) for severity in nature (based on mean annual temperature), hyphal growth, sporulation (weakly) and spore germination of the robusta population were all clustered around 23 °C, and their thermal response were roughly bell-shaped around T opt . In contrast, while the hyphal growth rate of the arabica population displayed a similar thermal performance profile and T opt to robusta (which also corresponds to previously measured T opt for growth in other Fusarium species ranging from 16 to 27 °C 21,42 ), all other measures displayed higher relative performance for arabica isolates at lower temperatures. In climate models, the arabica population displayed a lower inferred T opt and higher relative severity at lower temperatures. In the growth assays, sporulation and germination of arabica isolates occurred at a higher rate at low temperatures (10 °C and 15 °C) than robusta isolates. These findings support the hypothesis that the arabica population has higher tolerance to cold conditions observed at higher elevations where arabica coffee is grown 17,23 .
Our results imply that the arabica population is 'cold adapted'; which is observed in spore production and germination, which show more distinctive thermal responses. Potential explanations for this include accumulation of cryoprotectants 43,44 , alteration in composition of membrane lipids and fatty acids 44 , and evolution of psychrophilic (cold-active) enzymes in sporulating hyphae and spores 45,46 . Gene ontology terms for synthesis of trehalose, a commonly found cryoprotectant in microbes, have been found in the arabica population of F. xylarioides 47 , which might support the first explanation. As F. xylarioides is a soil-borne pathogen 5 , trehalose could prevent the arabica population from freeze-thaw damage which it potentially experiences in soils of Ethiopian Highlands during the cooler season 23 . Indeed, the arabica isolates had higher growth rates below 17 °C than the robusta isolates, a trend which reversed at warmer temperatures. Hyphal growth rates partly represent the infective abilities of the pathogen under ideal conditions, and could thus indicate a trade-off in virulence for increased resilience in cold environments. On the other hand, Coffea arabica is less resilient to diseases than robusta coffee 9 . Thus, arabica strains might need to invest less resources in hyphal growth to infect its target host plant compared to robusta strains. Both explanations suggest differential allocation of resources in growth and reproduction between the two populations, which could be investigated in future.
The phenotypic differences observed in the laboratory suggest that higher severity of coffee wilt disease (CWD) of the robusta population in Uganda and the Democratic Republic of Congo, compared to the arabica population in Ethiopia 17 , is not solely driven by different environmental temperatures across the countries, but also by differential thermal responses in growth and transmission of the two populations. However, growth, sporulation and germination cannot fully represent the life cycle of F. xylarioides. Other processes including spore dispersion and colonisation in host coffee plants will also affect its rate of transmission and infection. Therefore, it would be beneficial for future research to investigate the effects of temperature on these processes, especially on infection and growth rates in coffee plants, something we were unable to do. As F. xylarioides is soil-borne, future work could model investigate the effects of soil temperature (currently poorly measured at regional scales in Africa 48 ).
Through climatic modelling, we have also shown that although F. xylarioides is currently restricted to areas around the Congo Basin and Ethiopian Highlands, temperature patterns in several other coffee-growing regions of Africa are also suitable for one or both populations of the fungus. Many countries of these regions are important producers of arabica and robusta coffee, some relying heavily on coffee as a source of foreign currency 10 . These regions, such as west and south Africa (especially Angola, Zimbabwe, South Africa, Cameroon, Côte d'Ivoire, Gabon and Guinea), and Madagascar, could implement a range of measures, such as phytosanitary certificates; training extension workers to connect with smallholder growers; and conducting CWD surveys, to prevent future spread of the disease. These measures would ensure that any future outbreak is rapidly recognised www.nature.com/scientificreports/ and addressed. Several measures were proven highly effective in controlling the pathogen in Ethiopia and Tanzania, such as uprooting and burning infected coffee plants, phytosanitary certificate control, and education for coffee farmers on relevant information 13,49 , which can be promoted in other coffee-farming regions. Our future climate models predict overall decreases in severity because the mean temperature of modelled regions is predicted to be 27-35 °C by 2060-2080 29 , which is several degrees higher than the estimated optimum temperatures of both populations. Some countries are predicted to have climates favourable to increased severity, especially around the Congo Basin and Ethiopia. In principle, longer-term selection experiments could test whether the fungi are able to adapt to increase in temperatures, which then would alter predictions as the effects of particular temperatures on severity would then also change in the future. Though we have shown the strong effect of temperature on growth and transmission of both populations of F. xylarioides, temperature is not the only determining factor for CWD infections. Earlier research in other Fusarium species suggested the production and maturation of spores have a narrow range of humidity requirements 22 . Prior field observations have also shown alternations in rainfall promotes transmission of CWD, which provides the favoured humidity conditions for production of asci, and facilitate the release of ascospores from them 50 . Expanding our models to include other bioclimatic variables led to similar average severity predictions, but the regions projected to be suitable for F. xylarioides populations and subjected to future increases in severity were both more concentrated around equatorial Congo Basin (Fig. S2). However, both binomial generalised linear models had low coefficient of determination, especially for the arabica population (Supplementary Information S1), suggesting the need for more data and modelling to explain variation in their severity. Management approaches also strongly affect transmission of CWD. Mulching, herbicides, copper paints and hand-weeding suppress infection, while slashing weeds produces wounds on the plant, which allows more infection of F. xylarioides. Any systematic differences between regions or coffee species might further explain differences in severity between regions and host populations and would need taking into account. Thus, although we focused on temperature specifically, future research into how other climatic and agronomic variables affect the ecology of F. xylarioides populations and their corresponding coffee hosts would be helpful to predict future risk of disease across coffee-producing regions.
By demonstrating a match between macroecological patterns of disease severity and the thermal profiles growth and reproduction of disease-causing organisms in the laboratory, we provide evidence for temperature adaptation of CWD populations adapted to arabica and robusta coffee. As many parts of the world are predicted to experience warming, a combination of environmental modelling and physiological measurements with culture collection isolates can be a useful general tool to help to predict temperature-dependent life stages. Such information can feed into epidemiological models of future risk of plant diseases in changing climate.